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FOREWORD 


The information presented in this document was obtained as a result of 
the conical isogrid structural testing done by General Dynamics Convair 
Division under Contract NASA 8-29839. The contract was administered 
under the direction of John Key, Marshall Space Flight Center, NASA. 
Additional comparative data are presented from the cylindrical isogrid 
structural testing accomplished by NASA, Marshall Space Flight Center, 
in 1973. 

These structures were developed and fabricated by Convair under Inde- 
oendent Research and Development funding. Appreciation is expressed 
to Mr. Jack Furman of the Marshall Space Flight Center, NASA, for his 
contribution and continued interest in die development and application of 


isogrid to the aerospace structures. 
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APPENDIX C 


CYLINDRICAL ISOGRID STRUCTURES 
C.l INTRODUCTION 

Continuing work started on the conical isogrid adapter (see Appendix B) , the 
purpose of this handbook section is to present analysis techniques and data for eval- 
uating the load carrying capabilities of cylindrical isogrid structures subjected to 
compression and body bending loads. The results in this section are based on a 
structural test performed at MSFC on a Convair Aerospac3 furnished 120-inch diam- 
eter by 37-inch long cylindrical adapter. The loading fixture for this test was also 
designed and fabricated at Convair Aerospace. 

Figure C-la is a photograph and Figuie C-lb is a drawing of the flanged iso- 
grid adapter structure considered. Its construction is similar to that of the conical 
isogrid adapter except for the grid configuration which includes parallel-longitudinal 
rather than radial members. 

The approach of this handbook section will be: (1) to summarize test data 
resulting from a test on a cylindrical isogrid adapter, (2) to outline modifications to 
the techniques for predicting local and general instability described in Section B.2 
of Appendix B, (3) to develop a method for predicting extended local instability on the 
basis of beam column effects in two tandem longitudinal grid members , (4) to apply 
the extended instability analysis to predict general instability in the test specimen, 

(5) to demonstrate the effects of basic parameter and constant variations in the ex- 
tended local stability analysis, and (6) to evaluate effective skin width and compres- 
sion buckling of skin panels on the basis of test data. 

C. 2 SUMMARY OF TEST DATA 

This section contains a partial summary of the cylindrical isogrid structure 
test data for the final failure condition. Detail test results for all test conditions as 
well as a complete description of test setup and instrumentation is available in 
Reference C-l. 
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Figure C-lb. Isogrid Cylindrical Adapter 
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C.2.1 TEST SETUP AND INSTRUMENTATION 

Figures C-2 to C-4 are views of the test setup and Figure C-5 identifies 
the configurations and locations of the instrumentation used on the adapter and loading 
fixture. It is noted from these figures that there are 24 hydraulic loading cylinders 
and load cells to apply and measure loads on upper and lower loading fixtures assembled 
on the test specimen. The specimen is instrumented with about 370 strain gages and 
160 linear transducers. The strain gages measure grid member and skin stresses in 
the isogrid and flange transition areas of the adapter. The linear transducers mea- 
sure radial deflections at selected nodes. In addition to the above, eight linear trans- 
ducers measure movements of the upper loading fixture caused by load applications. 

The lower loading fixture is bolted to a rigid ground plane as is linear transducer 
support bracketry. 

As seen in Figure C-3 the linear ti'ansducers for measuring axial displace- 
ments are located close to the bolted flange between the upper loading fixture and test 
specimen. All linear transducers were kept in place up to and including specimen 
failure. Maximum transducer error is ±0. 00035 inch. 

The 24 loading cylinders act in 12 clevis assemblies as seen in Figure 
C-4. When a bending moment is to be applied the loading cylinders are programmed 
to apply loads that are nearly proportional to their distance from the axis about which 
the moment acts. This, in addition to an order of magnitude greater stiffness of 
loading fixtures compared to the stiffness of test specimen, causes the applied load to 
approximate perfect body bending. 

As seen in Figure C-5 instrumentation is concentrated in the vacinitv of 
the 0 circumferential station. This is the area where load peaking is generated and 
where initial specimen failure was expected. 

There is one jack bolt, as per Figure C-6., in each transition at the upper 
and lower bolted flanges (Figure C-6). Nut3 on these bolts are snugged up against 
the flanges to create direct load paths through the bolts and to minimize effects of 
flange bending or prying. Measured bolt tensile strength is 2100 lbs. 
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C.2.2 SPECIMEN FAILURE CONDITION 

Figure C-7 shows the test specimen carrying 250% of the maximal load 
for the failure loading condition 9 (see Ref. C-l for cylinder and load cell loads for 
this condition). Loading condition 9 includes combined axial and body bending loads 
to produce 636 lbs/in. specimen edge load intensity at circumferential station 0 for 
100% nominal load. For the 250% nominal loading or 1590 lbs/in. edge intensity the 
wrinkles in the skin are on ;he order of 0. 07-inch. 

The specimen f. tiled at close to 260% nominal load or 1654 lbs /in. Figure 
C-8 shows the failed specimen after it was removed from the loading fixture and Fig- 
ure C-9 shows a close up of a typically failed area. Large deformation damage in this 
photograph is attributable to post buckling effects. 

Initial failure probably occurred between circumferential Stations 54 and 
7. Damage beyond these stations is due to sudden load relief in the initially failed area 
and inability of the loading system to dump applied loads quickly enough to avoid secon- 
dary damage. 

Specimen failure was probably first sensed by a sudden loud report it 
created. This triggered a reaction in a matter of seconds by test personnel to dump 
load; however, this was not fast enough to avoid considerable secondary damage. 

C.2.3 INSPECTION OF DAMAGED SPECIMEN 

Figure C-10 diagrams damaged areas of the test specimen to indicate 
failed grid members and directions of structural deformations near damaged areas. 
Figure C-10 also lists measured or estimated grid- member and skin cross sectional 
dimensions in damaged areas. 

C.2.4 STRESSES IN FAILED GRID MEMBERS 

Figures C-lla to C-lld are plots of actual or representative inner flange 
stresses vs. precent load in some of the longitudinal failed grid members indicated 
on Figure C-10. Presented plots are based on availability of inner flange strain gage 
data. Since anticipation of all failed members was impossible, strain gage data is 
not available for many failed members. Strain gage numbers and corresponding 
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Figure C-lld. Stresses in failed grid members (longitudinal) 
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circumferential station numbers (per Figure C-5) are labeled with an asterisk in 
Figures C-lOa to C-lOd to indicate that stresses are for members representative of, 
but not the actual, failed grid member. When available these representative strain 
gages are chosen to have positions with respect to the upper bolted flange that are 
symmetrical with positions of failed grid with respect to the lower bolted flange. The 
tabulated cross sectional properties of failed grid members are average values based 
on the labeled representative grid members. 

Figure C-12 shows plots of available measured stresses vs. percent load 
in failed diagonal grid members. 

Figures C-ll and C-12 indicate different degrees of stress peaking or 
relief with load in the last 10 percent of load increment. From these plots tne most 
likely candidate for having precipated general instability is the failed grid member at 
circumferential Station 0. Apparently as seen in Figure C-llc this member is ex- 
periencing severe inelastic deformation in the last 10 percent load increment (i. e. , 
stress exceeds the 61,000 psi yield stress of 2024-T851 aluminum). This is not a 
certain cause of general instability since there are other measurements of comparable 
and considerably greater plastic deformation in areas showing no fractures. 


C.2.5 


STRESSES NEAR PEAK LOAD INTENSITY AREA 


Axial load (F) and bending moment (M) in a grid member cross section 
is computed from simultaneous solution of the following equations: 


and 
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where 

S , S 2 are respectively measured stresses on the inner flange and 
outer flange on skin side, 

C distance from neutral axis of grid member I-beam cross 

o 

section to extreme fiber on inner flange, 
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Figure C-12. Stresses in failed grid members (diagonal) 




s distance between extreme fibers on inner and outer flanges 

I, A section moment of inertia and cross sectional area based on a 

skin width equal to 47 times the skin thickness being effective 
in contributing to these properties. 

Based on average grid member cross section dimensions in Figure C-lb and the 

2 

computerized analysis described in Subsection C.3.2.3, C q = 0.504 inch, A = 0.162 in , 
and I = 0.013 in. 4 (s = 0. 73 inch). 

Using the above and available S and S measured (strain gage) data in Eqs. 2-1 

and 2-2 grid member axial forces and bending moments are computed for 200, 220, 

240, and 250% load conditions for part of the structure near the peak load intensity 

area. These are shown in Figure C-13. When measured stresses S or S exceed 

X 6 

61,000 psi (material yield stress) actual stresses are assumed to be 61,000 psi. Data 
in Figure C-13 are labeled with an asterisk to indicate measured stresses that exceed 
61,000 psi. Negative moments imply compression in outer flange. 

The following observations are made after an examination of data in Figure 
C-13. 

1. When the structure is primarily in its elastic range (i. e . , up to 200% 
load) axial forces in the diagonal members are about 10% of axial forces 
in longitudinal members. At 250% load diagonal members experience 
as much as 25% of the load in longitudinal members. This is caused by 
load redistributions due to local yielding in the structure. 

2. Positive moments in transition structure grid members at the upper 

and lower bolted flanges indicates possible failure of jack bolts in flanges. 

This failure would cause load to be relieved from the skin side of the 
structure and to be picked up on the inboard side of the loading fixture flange. 

The high bending moment in the failed longitudinal grid members at 
circumferential Stations 0 and 59 may have precipitated overall failure 
by failing inner flanges respectively in compression and tension. The 
negative bending moment in the failed grid member at Station 59 Is due 
to beam column effects while at Station 0 it is due to flange kick moments. 
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. Shear flow in skins and load&in diagonal members accounts for variations 
in loads along longitudinal sobers. 

4. Since diagonal grid members experience less load than in-line grid members 
it may be reasonable to assume that 70 effective skin width should be 
used in computing the cross sectional properties of diagonal members. If 
this were done the loads in the diagonal members would be greater than 
those indicated on Figure C-13 (by an estimated 17%). 

C.2.6 APPLIED LOADS 

It is assumed that loading fixture rigidity is an order of magnitude greater 
than that of the specimen. The 260%-of-loading-conditton-9-faihire load then pro- 
duces theoretical edge load intensi titles, N, and stresses in the transition structure. 


S^, at stations in the vicinity of expected peak load intensities. 
C-l was computed using: 

P 


N = 


21TR 


cos0 

trR 2 


See Table C-l. Table 


(2-3) 


where 

P is the axial load, 

M. body bending moment, 
b 

R cylinder radius, 

6 angle between Station 0 and reference station. 


Table C-l. Applied edge loads 


STATION 

N, lbs/in. 

S t , psi 

0 

1655 

22,761 

1, 59 

1650 

22,692 

2, 48 

1637 

22,514 

3, 57 

1616 

22,197 

4, 56 

1583 

21,771 

5, 55 

1544 

21,235 

6, 54 

1497 

20,588 

7. 53 

1442 

19.832 
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Measured specimen stresses In goemetrically identical lower-flange-transl- 
tion members which are in line with longitudinal members are plotted in Figure C-14 
as a function of percent loading. Station positions in Figure C-14 are projected on a 
line normal to the bending axis. The stresses, S^, from Table C-l , are also plotted 
on Figure C-14. The substantial difference between theoretical and measured stresses 
(at 260 and 250%) is due to toe following collective or individual cases: 

1. Reactions in the basic specimen structure produce bending moments in 
the transition structure. T1 3 measured stresses are the result of axial 
loads and flange bending moments. 

2. Crippling of jack bolts in loading fixture transfers loads to the inner flanges 
and creates a kick moment. 

3. Insufficient body rigidity in the loading fixture causes hydraulic cylinder 
loads to induce local load peaking. 

4. There is a discrepancy between loads that are proporational to distance 
from loading fixture bending axis and actual loads used in test. In both 
cases at 100% load, the total applied axial load is 120,000 lbs and body bending 
moment is 3.6 x 10® in. -lb. Figure C-15 shows differences between propor- 
tional and actually used forces for 260% applied load. The 883 lbs force for 
example is the difference between forces actually used and that would have 
resulted from proportional loading in hydraulic cylinders 1, 12, and 11. The 
force discrepancies are 4% for toe 883 lb force, 3% for 464 lbs, and 27% for 
2558 lbs. Directions of discrepancies are such that if the loading fixture 
rigidity were negligible compared to the specimen load peaking at the high 
load intensity area would be 15. 6% greater than predictable on the basis 

of a perfectly rigid loading fixture; it is estimated that the applied load 
error is not over 3 to 4% due to this cause. 

Effects of above factors influence the out of plane character of measured 
displacements in the upper loading fixture as shown in Figure C-16. Changes in 
fixture plane slope (in Figure C-16) is also influenced by local changes in adapter 
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Figure C-15. Differential forces on loading fixture 
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structure axial rigidity due to skin buckling. Skin buckling occurs at a fraction of the 
load carrying capability of the basic structure. This buckling is estimated to bo 
limited to the central 2/3 of the adapter length. 

The load peaking at circumferential Stations 0 and 56 as evidenced in 
Figure C-14 indicates the likelihood of failure having been precipated at these stations. 
General instability is most likely to have occurred at Station 56 where load peaking 
is greatest. 

The back-to-back strain gage data available on the longitudinal grid mem- 
bers between circumferential Stations 0 and 59 are used to compute applied axial (N) 
and kick moment (M) edge loading. These quantities are plotted vs. percent load in 
Figure C-17.The following assumptions were made to arrive at the computed values. 

1. The skin in the transition structure is fully effective. This results in 

42-times-skin-thickness effective skin width and resulting grid section 

2 4 

properties are C = 0.494 in, A = 0. 154 in and 1=0. 0124 in . 
o 

2. Axial forces and bending moments in grid members were computed per 
procedure described in Subsection C.2.5 and averaged for the four grid 
members involved. This assumes that the transition structure if fully 
effective in smearing loads. The highest peak edge intensity in the trans- 
ition structure at 260% load is 2224 lbs/inch, and the largest kick moment 
is 427 in-lb/inch. 

3. Two plots are shown for both N and M in Figure C-17. The lower values 
derive from assuming that the actually transmitted stress never exceed 
the 61,000 psi yield stress even though the measured gage streises are 
greater than 61,000 psi. The higher values are based on the measured 
stresses being the actual stresses even though they are over 61,000 psi. 

Data in Figure C-17 reasonably confirms the estimated 1654 lbs /inch 
failure load, however, this load is not necessarily at the location for which Fig- 
ure C-17 applies. 
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C.2.7 TOPOGRAPHICAL PLOTS 

Topographical and differential topographical plots of radial deflection 
contours in the test specimen, as shown respectively in Figures C-18 and C-19, are 
obtained through automated data reduction and plotting. Linear transducers measur- 
ing radial deflections at structural nodes (see Figure C-5) provide basic input for 
generating the topographical plots. 

C. 2.7.1 Program description 

Data reduction for isogrid topographical plots is broken down into four 
separate programs. Each program performs a unique function in the total process of 
conversion to final plot. 

The first program takes the GE 235 computer data tape produced by the 
MSFC-Structural Test and Checkout System (ST ACS) and converts it to 60 bit words, 
acceptable to the CDC CYBER 70 computer. 

The second program takes the new data tape and decodes the 60 bit words 
into real values which are the same values as those in the GE 235 tape. 

Program three writes the load scans on a tape for the plot program. This 
program b; way of input parameters selects the load scans wanted, and bv use of 
linear interpolation computes the values for non-existant nodes from the measured 
node values. 

The final program plots the contours from the data tape constructed from 
program three. The program can check for as many as thirty different contour levels 
that show up as solid lines for positive deflections and dash lines for negative. This 
program generates either straight plots and/or first or second order differencing 
plots which are determined by the input parameters. The actual plots are produced 
by the SC 4020 plotter. 

The smaller the selected increments between contour levels the closer 
are the contour plots and the more pronounced are surface slope indications. 
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C.2.7.2 Preliminary evaluation 

Examination of the contour plots (Figures C-18 and C-19) in fie vicinity 
of circumferential Stations 0 and 56 indicates that surface grades in thes* areas are 
large compared to surrounding areas and that much of this contour irregularity 
developed in the 240 to 250% load increment. This tends to confirm the conclusion 
in Subsection C. 2. 6 Ji&t failure probably started near circumferential Station 0 or 
56 and most likely Station 56. 

C.3 ANALYSIS OF TEST RESULTS 

The purpose of this section is to use test results to: 

(1) determine modification to the knock down factors for general instability 
described in Appendix B t 

(2) develop an extended local instability analysis and evaluate parameters 
in this analysis to obtain agreement with test results, 

(3) e vatuate the effective skin width included in grid member cross section, 

(4) develop relationships for predicting local compression buckling in isogrid 
skin panels. 

C.3.1 MODIFICATIONS OF GENERAL INSTABILITY KNOCKDOWN FACTORS 

U sing computer programs based on the analyses in Section B. 2 of Appen- 
dix B, Figures C-20a to C-20d were produced to compare actual specimen perfor- 
mance with four different predicted specimen performances. Figures C-20a to C-20d 
include the following curves: 

Curve labeled 1 is for general instability (i. e. , it designates a plot of 
allowable edge load intensity N , vs. node-to-node spacing a. for 
general instability as a limiting condition. ) 

Curve labeled 2 is for skin buckling (i. e. , it indicates the edge intensity 
at which skin buckling starts. Since the test specimen includes I-beam 
grid members skin stability is not critical to overall structural load 
carrying capability. This is similar to a condition common in skin stringer 
designs. ) 
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Curves 3 and 4 (not shown in Figures C-20a to C-20d) are respectively 
for web crippling and flange buckling. Edge intensities that create web 
crippling or flange buckling are greater than chosen scale values there- 
fore they do not appear on Figures C-20a to C-20d. 

Curves labeled with cross marks are for edge intensities at which material 
yield strength or grid member columr buckling is the limiting factor. In 
Figures 20a and 20b the cross marked curves indicate material yield pro- 
perties while in Figures C-20c and C-20d they indicate column buckling. 

The four predicted performances are: 

M-1A (Figure C-20a) in which general instability is based on the McDonnell 
Douglas analysis (see Appendix B) and assumes that the skin is fully 
effective in reacting load (i. e. , even though the skins elastically buckle 
they develop tension fields that cause them to be fully effective in reacting 
loads). 

M-1B (Figure C-20b) in which general instability is based on the McDonnell 
Douglas analysis assuming that a width of skin equal to 47 times its thick- 
ness, t , effectively contributes to the cross sectional properties of the 
I-beam grid member (e. g. , 23.5^, on each side of outer flange). 

B-1A and B-1B (Figures C-20c and C-20d) which are respectively like M-1A 
and M-B except that the general instability curves in this case are based 
on a Convair Aerospace analysis presented in Subsection B. 2. 3 of Appen- 
dix B. 

The structural parameters used in the above analyses are based on aver- 
age values in Figure C-l and listed on Figures C-20a to C-20d by the H. P. Calculator 
Plotter (see Section B.2. 1 of Appendix B for definition of symbols). 

In all cases predicted structural performance is dictated by general in- 
stability for the 7. 338 inch node-to-node spacing in the test specimen. 
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It is clear from Figures C-20a and C-20c that the assumption of fully 
eriective skin causes an excessive discrepancy between predicted and actual structural 
performance. These analyses are assumed invalid in evaluating test data. 

From Figure C-20d it is seen that predicted performance based on gen- 
eral instability and column buckling are within about 5% of each other. This is within 
range of theoretical accuracy. Therefore the knock down factors corrections in 
Table C-2 are recommended to bring theory in line with test data. 

Table C-2. Summary of theoretical and actual critical loads 


Damaged area: 

Average actual failure load: 
Theoretical failure loads 
General instability 
MDAC * analysis : 
Convair** analysis: 
Column budding: 
Knock down factor correction 
General instability 
MDAC analysis: 
Convair analysis: 
Column buckling: 


Center of panel No. 1 
(see Figure C-lb) 

1655 lb/in 


2140 ib/in 
1770 lb/in 
1854 lb/in 

theoretical/N actual) 
CK 

1.29 

1.07 

1.12 



*per Subsection 8.2. 2 of Appendix B 
**per Subsection B.2. 3 of Appendix B 

Insufficient data is available at this time to preclude whether general 
instability or column buckling is the dominant failure mode. 
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C.3.2 EXTENDED LOCAL INSTABILITY ANALYSIS 

Beam column instability of two tandem longitudinal grid member* will be 
considered as a basis for evaluating extended local load carrying capability of iso- 
grid structures. Grid orientation "b" in Figure B-4 of Appendix B is assumed for 
purposes of this analysis. 

C.3.2.1 Description of model 

The model to be analyzed includes two initially straight members the 
neutral axes of which intersect at a point offset e with respect to a line along which 
applied forces F act as shown in Figure C-21a. The offset is taken as positive when 
it is in the sldn-to-inner flange direction. It can represent anticlastic or synclastic 
curvature, misalignment of applied forces, and eccentricity between applied forces 
and beam neutral axis. 

A fixity due to a torsional spring rate is assumed to exist at points 1 
and 3 in Figure C-21a. This spring rate, about a normal to grid-member -plane 
axis, is determined from the model in Figure C-21b which includes an in-line grid 
member of arbitrary length a". Likewise die linear spring restraint k at point 2 
(also used in the model for determining k$) is based on the model in Figure C-21c. 

An arbitrary grid member length a' is used in this model. The equivalent lengths a' 
and a" are obtained from: 

a’ = K f a (3-la) 

a" = K & a (3-lb) 

where factors K, and K are determined to equate theoretical and measured data. 

I a 

Isogrid grid member cross section shown in Figure C-21d is taken normal 
to the neutral axis and to include 47 effective skin width contributing to the cross 
sectional moment of inertia about the indicated neutral axis. 

While the analysis to be presented is of a flanged isogrid structure it is 
readily extended to: 

(1) Unflanged isogrid with skin, 

(2) Flanged isogrid wit' <xit skin, 
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Figure C-21,Beam column mod* i for extended local instability analysis 
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(3) Frame stabilized skin-stringer structures. The choice cf K and K 

ft X 

factors would have to account for differences in frame and stringer cross- 
sectional properties if any, 

(4) Flanged or unflanged waffle structure. 

C. 3.2.2 Definition of symbols 

The symbols used are as given in Subsection B.2.1 of Appendix B, those 

above, in Figure C-21 and the following: 

= allowable edge load for conditions that yield stress F^, 

is not exceeded In the structural model in Figure C-21, 

A = A' = area of equivalent I-beam member cross section Including 
47 1, width of effective skin. 

So, Sa = grid member stresses at x - 0 and x = a in F. re C-21a, 

I = is about the neutral axis shown in Figure C-21d, 

b = (b x + b 2 )/2, 

Tj = the fraction of the total losd reacted by in line grid members. 

C.3.2.3 Isogrid cross sectional properties 

The properties A, I and Co are computed from the following on the basis 
of the parameters shown in Figure C-21d. 

A = (t ; + c) w + 47 1, 2 + (s - c + tj b (3-2a) 

a 1 2 

t ^ 3 . ,2 , 3 X 

1 = 12 w c c w (C ^ 2) + b (s - c - y Yz 

+ ( s-c-t 2 )b[c + i(s-c-t 2 )] Wt 2 

+ t 2 w (s - t 2 /2) 2 + Y 2 < 47 V ‘i + 47 \ < 8 - \ /2)1 (3— 2t>) 

Co £wc - | (2c) + 47 1^ + wt ? + sb - t g bj - ~- 

+ 47 t * (S - tj/2) + wt 2 (S - t 2 /2) (s - t 2 ) 2 (3-*. ’) 
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C. 3.2.4 Analytical Formulation 

Aaanmtng a straight line grid mem 'jer between points 1 and 2 in Figure 
C-21a and that y Is a deflecticn with respect to ‘dds line 

A 


M = -El 


da? 


(3-3) 


where M, and applied bending moment, is positive when it muses compression cm the 
right side of the grid member. 


Between paints 1 and 2: 


-k e -y 
a t 1 J 2 2 


(3-4) 


where 0^ is the an gular deflection at point 1 with respect to die line between points 


1 and 2 and is positive when clockwise relative to this line. y g is die displacement at 


point 2. 

From Eqs. 3-3 and 3-4 


El 


$--*4? -T 1 ]■•«.*, 


(3-4a) 


The following is a model far Eq. 3-4a: 


y = ^ y + x + ^ 


(3-4b) 


Using La Place Transform techniques the following solution has been developed for 
Eq. 3 -4b. 


y= - 


e. 


(-bj) 


1.5 


-(-bj) 


0. 5 


sin (x ) 


a i X fl 

b + b 
1 1 


cos (x y=B 1 ) -1 


(3-4c) 
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This is also a solution to Eq. 3 -4a when the following substitutions are made 


b - k y 
1 11 y 2 


k 14 e i 


1 -Fe 
b - ■ 

1 aEI 


= -k/2EI 


‘m'V* 


Eq. 3-4c is valid when the force, F, is compressive as shown in Figure C-21a and 
therefore b^ in Eq. 3-4d is negative (i.e., - b, is a positive value). 

Slope, e , at any point on dm elastic curve as a function of x is given by 


-v±/ 


(3-5) 


Using Eqs. 3-4 and 3-4d in Eq. 3-5 and performing integration 


6 = 0 


a x e x [cos (x y^)-lj 


+ _I_ J) J_L_\ FZlll 

EI ( (-b^) 1 * 5 (-b/* 5 ) (-b/* 


l l X k 14 ®1 *** (X k 14 G l x 


b l <-v 


_ ky n 1 2 

-Fe _2 x . ^ 

+ + — +k © x 

a 2 2 t 1 


Making use of the boundary condition 0j= 0 at x = a in Eq. 3-5a, ©2 is determined 


as follows: 


a - 1 .1 . F ea 

6 1 _ k EI ” b l k 7 2 ‘ y 2 

« 


(t - - k n n) 
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where 


V*‘ 


F jcos (a 7-bj) -lj 


Elb 


k a k aF 

-L. + _I£_ 

EI b El 


Fk^ sin (a/^bj) 


EI b x (-bj) 


0.5 


k 7 = 


-F 


cos (a^) -1 a 2 


2 


2 b. 


Making use of the boundry condition y = y at x = a and 9 from Eq. 

A X 

3-4c: 

b 1 k + k k /k EI 
1 71 72 74 2 


y 2 = 1-k 


75 


where 


k 


sin (a/^) 


71 


<“V 


1.5 


72 


.1 i. Fea 

-b k + 

17 2 


74 


sin(a/^b) k 

— * — |coa (a/^) -l] 


(-V 


ka 

k = ~~ - k k 
73 4 11 7 


. k k 73 k 74 

k 75 ‘' k ll k 7l" k 2 EI 
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3-6 in Eq. 
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(3- 7a) 

(3- 7b) 
(3- 7c) 

(3-7d) 

(3-7e) 
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Deflection y is numerically evahiatable from Eq. 3-7. y is theu used in Eq. 
3-6 to determine 6 . 

When the applied force is in tension, and therefore b in Eq. 3-4d is positive, 
the sine and cosine terms in die above expressions are replaced by hyperbolic func- 
tions which are expressed as: 


cos (a ^/-bj) = 1 1 exp (a/t^) +exp (-a y/bj j , 


and 


sin (a yfibj = ~~ [ exp (a JbJ - exp (-a/^) ] . 


Imaginar y terms resulting from positive cancel out when these substitutions 


are made. 


The equivalent spring restraints k and k in Figure C-21b and C-21c are 

* 


k = 


24 El 

(a'F 


(3-8) 


and 


l k/3 + EI/(a") I 


(3-9) 


where a' and a" are defined in Eqs. 3-1 and 3-2. 

When y and 6 are determined M can be solve-' from Eq. 3-4 and beam stresses 

2 X 

are found from 


MC 


S = 


M(s-C o ) 

° r S " *A 


(3-10) 


Maximum stresses occur at x - O, S 0 , and x » a, S a . 

The analytical approach is to determine the force F at which either S or S is 

o a 

at the yield stress F . Then 
cy 


N 7 -jos 60°) V a 


(3-11) 


t 
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C. 3.2.5 Computerized Analysts 

The analysis outlined in 9 bsection C.3.2.4 has been programmed for 
the Hewlett Packard 9810 Desk Computer and 9862A Calculator Plotter. Figure C-22 
is a flow diagram for this program. It has the following features: 

1. The inputs include parameters needed in the plot program as well as in 
Hie structural analysis part of the program. The latter includes an initial 
applied force, F-start, and an arbitrary increment 4 F. 

2. Initially A, I and Co are computed. The quantity is then computed and 
its sign is checked by an IF statement which decides whether hyperbolic 
or trigonometric functions are applicable in the analysis. Results of 
analysis using either of these functions are branched to LBL 1. 

3. Quantities y and 6 are next computed after all the k factors needed 

m X 

in their evaluation have been determined. 

4. It has been found that near critical loads a small increment in F causes 
radical increases in y 2 and 0^. If F is increased beyond the critical 
load point y quickly rises to positive and then negative infinity and be- 
haves like a tangent function for argument at */2 radians. To avoid 
resulting irrational solutions, especially when F start is chosen to be 
larger than experience with a few runs would indicate it should be, a test 
is performed to determine if y is greater than 0. 01 a. The 0. 01 is 

A 

arbitrary but found adequate to compute forces F within 1% of critical 
values. When the y > 1 0. Ola | statement is satisfied N is computed 

t* ( 

and sent into storage for plotting. Computed N^, y^. So, Sa and F are 
printed out at this time. 

5. If the y > 1 0. Ola \ test is not satisfied So and Sa at the two ends of the beam- 
column are computed and individually chocked against Fey. If Fey is not 
exceeded F is Incremented by AF and ‘he loop beyond LBL-?- is recycled. 

When Fey is exceeded at either end of the beam column N ? is computed and 
directed to LBL X 2 beyond which the indicated parameters are computed 

and printed. 
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Figure C-22 Computer program flow diagram 
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6. After a is used beyond LBL X 2 it is incremented by 1 in. and returned to 
LBL-r for program recycling. The program ends when a = a min. at which 
time tiie plot program generates vs. a. 

7. Force, F, is incremented from small to large values and the initial node-to-node 
spacing, a, is the largest to be plotted. This causes the m a ximu m F for any 

a to always a good initial value for the succeeding a (since the load reacting 
capability of a column either increases or remains unchanged as its length 
decreases). 

C. 3.2.6 Numerical Evaluation 

The following program parameter inputs include tj» Kg.> and Kj values 
found to produce analytical results that agree with the test results. 


r\ 


0.9 

w 

— 

0.415 

e 

= 

0. 02 in.. 

c 

= 

0.082 

s 

= 

0. 73 in. 

K 

a 


1.4 

AF 

= 

100 lb. 

k » 


1.4 

F-start 

= 

2000 lb. 

b 


0.053 

F 

rr 

60, 000 psi 

t. 

= 

0.040 

cy 


7 

1 



E 

= 

10 psi 

*2 

— 

0.052 


Effective skin = 47 

The V = 0. 9 was estimated from data in Subsection C. 2. 5. The e = 0. 02 
in. was computed on the basis of 0. 15 in. total anticlastic curvature offset measured 
in the foil length (37. 25 in. ) of the test specimen. The s, w, c, b, t and t values 

1 M 

are for the average grid member cross section from Figure C-lb. F- start and AF 

are arbitrary program operational values. F and E are for the 2024-T851 aluminum 

cy 

from which the specimen is fabricated. Finally K and K factors were determined 

& I 

by exercising the program. 

Influences of the different parameters affecting computed and automated 
plots of N ? vs. a are studies in Figures C-23a to C-23d. In these figures program 
input parameters are as listed above except for one parameter being studied and 
indicated on each figure. 
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Figure C-23d Comp 
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Figure C-23a shows the effect of changes in on N^. Note that the 

vs. a curve passes through the 1655 Ib/in. test value when Kf = 1.4. Figure C-23b 

shows the effects of changes in K . It is apparent that K is more significant than K . 

at a 

Note also that the N for K = 0. 1 is slightly less than it is for the less-end-fixity 
7 a 

case = 100. This is because in the first case peak stresses develop at x = o where 
the distance to the extreme fiber Co is greater than s - Co for the latter case in which 
peak stresses are at x = a. 

For values of a greater than 9 inch N ? values in Figure C-23b are con- 
trolled by the y >|0. 01a| test in Figure C-22. 

Figures C-23c and C-23d show effects of variations in e and effective 
skin. It is seen that e is not as critical a factor as the effective skin based on values 
chosen for each. 

Table C-3 is a computer printout for the inputs listed at the be ginning of this 

subsection. It is noted from observation of y and a values that N , So, Sa and F are 

2 7 

determined from the y >| 0. Ola | test for values of a = 15 to 10 in. For these values 

the full stress F of the beam column is not developed. For values less than a = 9 in. 
cy 

the full yield stress is developed at one end of the beam column (and y <0. Ola). 

C.3.3 EFFECTIVE SKIN WIDTH 

The object of this subsection is to use isogrid adapter test data to develop 
relationships for predicting the effective skin width in compression to be included in 
grid member cross-section. 

Preliminary evaluation of isogrid test data indicated that effective com- 
pression skin width acting with the stiffener is higher thf . predicted by the standard 
method used for skin-strl'-ger design. To effectively use this added area, a method 
for predicting effective skin width consistent with available test data is developed. 

The simplified analysis most often used for skin- stringer design assumes 
a plate of width 2 w simply supported on all four sides working at the same stress 
level as the stringer. 
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Table C-C. Computer printouts. 


0 . 9080000 ^ 
6.0200000 
0.7300060 
100.0000000 
2.000000000 63 
6.000060600 04 

*1.000000000 07 ^ TM13TTT®1 
0.4150000 f* INPUTS 
0.0320000 
1 . 4000000 
1.4000000 
0.0530060 
0.0400900 
0.0520900 J rn 



9. 1623930- 

0.0129647 — I 


397.9291763 
0. 1806909 
2. 417483474 04 
2.349621992 64 
3.60 8 0 0 0 0 0 0 0 3 
15.0080800 

3 f‘5 . 7 469 1 15 



0. 1645130 
2.753163533 04 
2.677172763 64 
4.100O08000 03 
14 . 0800600 


746.4949728 
0.1153340 
4.261709909 04 
4.154176636 04 
6.400000000 03 

11.0000000 

987.9394406 
0.1101646 
5.133463218 94 
5.004995550 04 
7.700000900 03 

10.0000000 

1 . 297294215 03 
0.0815181 
6.071331274 04 
5.929973891 04 
9.100000000 93 

9.0000000 

1.491531948 93 
0.0334215 
6.091579195 04 
5. 842004O83 84 
9.300000000 03 

8.0000000 


1.741266175 83 
0.O181030 
6.063564*649 04 
5.964494402 04 
9.500000000 03 

7.0000000 


4 6 3 • 8676694 
9.1452583 
3.147185460 04 
3.0628". ^79 04 
4.700000000 03 

13.0000000 


2.031477295 93 
0. C' 03138 
6.928889269 94 
5.937.137530 04 
9.590000008 03 
6.0088898 


583.O591903 
0.1411 165 
3.6951546 9 84 
3.596081 '"92 84 
5.500000000 03 
12. 000 CU00 


2.437772646 03 

0.00602:2 

6.099488248 94 
5.922411321 04 
9.500900008 83 
5.8009008 
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For the skin section in Figure C-24 critical buckling stress; 
/ t \ 2 

F. = KE (“ — ) 
c CR \2w / 

For 1 >>2w, K = 3. 62 for all sides simply supported. 

mu* 


r‘ JT 

i== t r r l ■ 
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Figure C-24. Typical skin stringer geometry 


Solving Eq. 3-12 for w in terms of F^ with K - 3. 62 


w 


- o-ntjf 


(3-12) 


(3~12a) 


This is the classical equation for effective skin width. For a typical 
wiuminum alloy working at F c = 40,000 psi, w = 15 t or 2w = 30 t is often used for 
preliminary sizing. 

C.3.3.1 Test Data 

One of the skin panels on the test specimen was instrumented with several 
rosette strain gages as shown in Figure C-25a. Figure C-25a and C-25b outline strain 
and stress data resulting from these gages. See Ref. C-l for detailed test data. 

Figure C-26 is a plot of data in Figure C-2Ca. Figure C~27 shows cal- 
culations of average stress and factor C in effective skin determination based on test 
data. From Figure C-27 the effective skin width is 


w ' 1 - 67t \/§'. 


(3-12b) 
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GAGES AT 1, 2 «t3 WERE USED TO 
CALCULATE THE STRESS 
DISTRIBUTION IN THE SKIN 
AND STIFFENER. 


( ) INDICATES GAGE AT 

RIGHT ANGLE TO LOADING 
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GAGES IN DIRECTION OF LOADING 
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-1776 
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-2969 
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3373 

-3033 

-3883 
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Figure C-25a Typical local specimen strain stress data 
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GAGES AT RIGHT ANGLE TO LOADING 


Gage 

Strain — p in. /in* (Cond 9A) | 

80% 

120% 

160% 

200% 

240% 

3091 

705 

1026 

1339 

1656 

1876 

3094 

-245 

-41 

209 
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70S 
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1499 

2772 

4123 

5821 
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721 

668 

475 

-599 


AVG STRAIN 


PT 




Strain — 

pin. An. 

(Cond 9A) 
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240T 


€ i 
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«• 
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mm 


a 

-1415* 
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1733* 

-5560* 

2399t 
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3255t 

2 

-1286 
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-23 Ou 
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-3301 

1720 

-4327 

2299 

-5467 

2611 

3 

-565 

230 

| -644 

492 

-746 

774 

-851 

1055 

-952 

1289 


Calculate stresses parallel and perpendicular to direction of loading* 

*■ + V€ X> - 77 ^ ‘*1 ^ I '*!) 


AVG STRESS 
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- psi (Cond 9A) 
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240*7 

PT 

*> 

a 


a 






°L 

1 1 

* 

-14503 

1181* 

-26960 

3451* 

-40155 

6152* 

-55849 

8438* 

-86648** 

8186* 

0 

-13017 

1617 

-22696 

4846 

-32134 

8419 

-41968 

11548 

-54042 

11202 

» 3 

-5723 

69$ 

-5727 
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-5923 

6348 

-6167 
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-6522 

11577 


-02 

* f H - € Skin ' 773 ^SWn “ c Stiff> 

* ■ Based on E - 10. o * 10®. Material has obviously exceeded the yield stress so this value is too high. 

•There were no rosette gages on the stiffener at point 1, so average stress is calculated using stresses 
perpendicular to loading direction at point 2. Assume stress perpendicular to ’oading direction at 
point 1 f(7j_) is equal to 
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Solving for « A : 

Frr^urp C-25b, 


tr A )( 


l-i/- 

L_ 


) - V( t 


Typical local specimen strain/ stress data. 
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DISTANCE FROM <£ STIFFENER-IN. 


Figure C-26 Grid-skin stress distribution 
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| AVG SKIN STRESS (KS1) AND EFFECTIVE WIDTH (IN) | 


80% 

120% 

160% 

200% 

(SAvg 

8.2013 

11. 1624 

14.0504 

16. 9290 

w 

1.681 

1.231 

1.040 

.901 

c 

1.64 

1.64 

1.69 

1.73 


Avg 




2.9725 a 


EFFECTIVE SKIN WIDTH, w = 


Avg 


a Stiff 


WHERE (5^ tiff * (pfl @ PT (1) FROM FIGURE C-25b. 



SOLVING FOR c 


t VeTF 

c 

WHERE t = .038 IN 
F c = tfStiff 

FROM THE ABOVE C Avg = 1.67 


Figure C-27 Effective skin claculations 
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C.3.3.2 Analytical Estimation of Effective Skin Width 

From the analysis of compression budding of skin panels in Subsection 
C.3.4 K c for the reference panel is 17.25. 

Kg ir 2 

K = — — - 15.59 (3-13) 

12 (1 - v 2 ) 



(3-14) 


and since K - 15. 59 


w = 1.97 t 



(3 -14a) 


This is obviously optimistic since the K is based on all edges having substa nti al support 
above simple support. For all edges simply supported K = 6.78 and 


w = l 


• 3ot \S 


(3-14b) 


Obviously the actual value of w falls somewhere between the values defined by equations 
3-14a and 3-14b. Therefore letting 



The effective value of c can arbitrarily be defined as 

c 0 - 1.30 

c = 1.30 + 

2 

and the effective skin width is 



Using Eqs. 3-14d and 3-14e for the test panel 


(3-14c) 


<3-14d) 


<3-14e) 



„ „ A 1.97 ■ 1.30 , Ai 

c - 1.30 + r = 1.64 

2 

which agrees closely with the measured values. 

C.3.3.3 Analysis Procedure 

1. Calculate skin panel buckling coefficient, K c , using procedure outlined in Subsec* 
tion C.3.4. 

2. Calculate K = 


Kc * 2 


12 a - V 2 ) 


3. Calculate c. 


JK 


4. Calculate c = 1.30 + 


Cq * 1.30 




5. Calculate effective skin width 
w = ct^ 

As a conservative approximation w can be defined as 

w - 1.30t./-£ • 

V F C 

Using this relationship for w and F c = 40,000 psi (aluminum) 

w = 21.06 t and 2w r 42. 12 t (compares to 2w = 30 t used for skin-stringer) 
which can be used for preliminary isogrid sizing. 

C.3.4 COMPRESSION BUCKLING OF SKIN PANELS 

The objective in this subsection is to use test data to develop relation- 
ships for predicting compression buckling of isogrid skins. 


C.3.4. 1 Measured Buckling Data Evaluation 

Buckling of skins was measured by back-to-back rosette strain gages 
on four triangular sKin panels. Outputs of these gages are plotted in Figures C ?8a 
to C-28d for loading condition 7 in Reference 1. The figures also showc and €, 
computed from this data. 

€ is the average strain in the direction of loading. € is the average 

strain perpendicular to the direction of loading. The point of incipient buckling is 

C-60 



I 


defined as the point at which the curves in Figures C-28a to C-28d start to deviate 
from a straight line. Critical buckling stress. 


a 

CR 


2 1 + ^ C 2^ 


(3-15) 


€ and £ being the measured strains at incipient budding 
X M 

From Reference C-2 critical buckling stress can be defined as 


<r 


CR 


sili (4 

12(l-/i 2 ) ' b ' 


(3-15a) 


solving for K (the budding coefficient) 

c 2 

12(1-// ) <r CR 

c ir 2 E (t/b) 2 

where 

t = average skin thickness, inch 
b = panel size, inch. 


(3-15b) 


K calculations are summarized in Table C-4 using data from Figures C-28a to 
c 

C-28d. From Table C-4 the average value of K is 17. 0. 

c 


Table C-4. K Calculations 

c * 


Panel 

Gages 

t 

(in) 

b 

(in) 

6^ Micro 
Strain 

€ Micro 
2 Strain 

a CR 

(psi) 

K 

c 

1 

3055/3058 

3057/3060 

0. 039 

6. 35 

-580 

180 

-6069 

16. 95 

2 

3091/3094 

3093/3096 

0.038 

6. 35 

-550 

240 

-5861 

17. 25 

3 

3121/3124 

3123/3126 

0.37 

6. 35 

-520 

160 

-5446 

16. 90 

4 

i 

1 

3147/3150 

0.037 

6. 35 

-520 

150 

-5480 

17.00 
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gure C-28b Skin strains vs load 
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Figure C-28c Skin strains vs load 
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Figure C-28d Skin strains vs load 
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C.3.4.2 Effects of Skin Edge Fixity 

Actually Kg is a function of the rigidity ratio, 

c = 2*22X1 kin ( 3 - 16) 

®max grid member 

where © max is the maximum edge rotation for a trian gular skin panel subjected 
to a unit distributed edge moment and © max member is the maximum rotation for 
a grid member subjected to a unit distributed moment. 

Using Table C-5 and Figure C-29, rigidity ratios are calculated in Table C-6. 

where 

I = Miminum Moment of Inertia about the 0-0 axis in Figure C-29, in£ 

4 

J = Polar Moment of Inertia, in. 


Table C-5. Panel Geometry Values 


Panel 

Cages 

Panel 

t 

d 

£ 

w 

c 

b i 

b 2 

»2 

1 

3055/3058 

3057/3060 

0. 039 

6. 35 

7. 33 

0.415 

0. 085 

0. 055 

0.083 

0.057 

2 

3091/3094 

3093/3096 

0.038 

6. 35 

7. 33 

0.415 

0. 085 

0.054 

0. 083 

0. 052 

3 

3121/3124 

3123/3126 

0. 037 

6. 35 

7. 33 

0.415 

0. 087 

0. 058 

0. 086 

0. 056 

4 

3147/3150 

3145/3148 

0.037 

6. 35 

7. 33 

0.415 

1 

0. 085 

• j 

0. 053 

! 

0.083 

0. 051 



Figure C-29. Average Panel Geometry Sketch 
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Table C-6. Panel Rigidity Ratios 


Panel 

Gagee 

.. 

D 

I 

J 

.. 

a 

Gmax plate 

3max stiff 

€ 

1 

3055/3058 

3057/3060 

57.04 

8.611x10“^ 

1.720x10*4 

1.323 

2,78x 10” 2 

6.91xl0* 4 

40.2 

2 

3091/3094 

3093/3096 

52.76 

8.313x10-4 

1.656x10-4 

1.325 

3.01x10-2 

7.18X10* 4 

41.9 

3 

3121/3124 

3123/3126 

48.70 

8.691x10-4 

1.847x10-4 

1.283 

3.26x10-2 

6.37X10* 4 

51.2 

4 

3147/3150 

3145/3148 

48.70 

i — 

1 8. 251 x IQ -4 

1.636X10" 4 

1.328 

3.26x10-2 

7.27X10* 4 

44.8 


gw 3 ,375 bj 3 dbg 3 t 2 w 5 
12 12 + 12 * 12 


wc 3 ,375bi® db'"* wt 2 3 


In Equation 3-16 e is computed witii the following. 
From Reference C-2 


max s idn ' 4D 


where 


D = 


Et* 3 


12(1"M 2 ) 


from Reterence C-3 


e 


max grid member 2 JG 


~ - a tanh 
4 


(A) 


where 


a = j/EI/JG 


h = Stiffener height 
jt = Stiffener Length, in. 


<3-16a) 
(3 -16b) 


(3-16c) 

(3-16d) 


Triangular panel compression buckling coefficient as a function of panel to edge mei. 
ber rotational rigidity ratio is plotted in Figure C-30 fr- m the data in Tables C-4 ano 
C-6. 
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C. 3.4.3 Analysis Procedure 

1. Calculate rigidity ratio, c as outlined. 

2. Using the calculated value of c obtain K(. from Figive C-30. 

3. Calculate critical buckling stress (uniaxial compression) for the panel from Equa- 
tion 3-lSa. 

C.** SUMMARY AND RECOMMENDATIONS 

; . A limited treatment of the test data has been presented. The data has been used 
to evaluate knockdown factors for general instability, extended local instability, 
effective s 1 ^ width and compression budding of skin panels. 

2. Failure of the loading fixture bolted flange jack bolts probably produced premature 
failure of the test specimen. It may also explain failure along the adapter lower 
end, althoug h such failure modes in cylinders are common. The load reacted by 
the cylinder would likely have been 20% higher, were it not for the flange failure. 

3. A unique automated method has been developed for generating topographical pilots 
of radial deflection contours in the cylindrical adapter. The plots are used to 
anticipate and identify areas of incipient and actual structural failure. 

4. A developed extended local instability analysis generally agrees with the general 
instability analyses presented in Appendix B. 

5. Effective skin width and skin budding evaluation techniques have been developed 
for isogrid structures. 

6. It is recommended that test data be used in further evaluations including discrete 
element modeling of the test specimen. 

The knowledge gained in conducting and evaluating the isogrid cylinder test will be 
significant in planning and executing future isogrid cylinder tests. It is recom- 
mended that if the subject test be rerun: (1> toe loading fixtures be stiffened and 
strengthened, to ensure two orders of magnitude greater stiffness and strength In 


f* 

ft* 
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the fixture as compared to that In the adapter; (2) greater concentration of back- 
to-back grid member strain gages be used In the expected failure area; and (3) 
only body bending load application be used, to maximize localization of failure 
area. 
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